Cosmic string induced sheet like baryon inhomogeneities at quark-hadron transition 
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Cosmic strings moving through matter produce wakes where density is higher than the back- 
ground density. We investigate the effects of such wakes occurring at the time of a first order 
quark-hadron transition in the early universe and show that they can lead to separation of quark- 
gluon plasma phase in the wake region, while the region outside the wake converts to the hadronic 
phase. Moving interfaces then trap large baryon densities in sheet like regions which can extend 
across the entire horizon. Typical separation between such sheets, at formation, is of the order of a 
km. Regions of baryon inhomogeneity of this nature, i.e. having a planar geometry, and separated 
by such large distance scales, appear to be well suited for the recent models of inhomogeneous nu- 
cleosynthesis to reconcile with the large baryon to photon ratio implied by the recent measurements 
of the cosmic microwave background power spectrum. 
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I. INTRODUCTION 

Generation of baryon inhomogeneities during quark- 
hadron transition is an important issue due to its con- 
sequences for the nucleosynthesis, and the possibility of 
creating compact baryon rich objects (lj. Recently, there 
has been a renewed interest in the investigation of nu- 
cleosynthesis in the presence of baryon number inhomo- 
geneities, the so called inhomogeneous big bang nucle- 
osynthesis (IBBN) [^|J|]. This originates from the re- 
cent measurements of the angular power spectrum of the 
cosmic microwave background radiation (CMBR) by the 
BOOMERANG, and MAXIMA experiments [§, which 
lead to a value of baryon density to critical density ratio 
f2f, which is larger than the value allowed by the conven- 
tional models of homogeneous nucleosynthesis. Earlier, 
there have been numerous investigations of the nature 
of baryon inhomogeneities generated during a first or- 
der quark-hadron phase transition pi and its effects on 
the abundances of various elements |6[. In these inves- 
tigations, baryon inhomogeneities arise due to moving 
bubble walls at the transition, with baryons getting con- 
centrated in the remaining localized quark-gluon plasma 
(QGP) regions. Typical separation, linhm-, between such 
localized baryonic lumps is of the order of separation be- 
tween the nucleation sites of the hadronic bubbles, which 
is at most of the order of few cm at the end of the quark- 
hadron transition for homogeneous nucleation ^|J^]. In 
order that these baryonic lumps survive various dissipa- 
tive processes, such as baryon diffusion, neutrino infla- 
tion etc., up to the time of nucleosynthesis, the required 
separation of these lumps is of the order of few meters 



There have been some attempts of getting a larger sep- 
aration between baryon inhomogeneity regions, linhm-, 
by deviating form the homogeneous nucleation theory. 
Christiansen and Madsen have argued that sufficiently 
large values of linhm can arise in the presence of impu- 
rities due to heterogeneous nucleation of hadronic bub- 
bles. It is mentioned in ref. that possible sources of 
such impurities could be primordial black holes, cosmic 
strings, magnetic monopoles, or relic fluctuations from 
the electroweak scale. Hadronic bubbles are expected to 
nucleate at these impurities with enhanced rates, lead- 
ing to separation between baryon rich regions being of 
the order of typical separation between impurities. Re- 
cently, Ignatius and Schwarz have proposed || that the 
presence of density fluctuations (which could arise at the 
time of inflation) at quark-hadron transition will lead to 
splitting of the region in hot and cold regions with cold 
regions converting to hadronic phase first. Baryons will 
then get trapped in the (initially) hotter regions. Esti- 
mates of sizes and separations of such density fluctuations 
were made in ref. jq] using COBE measurements of the 
temperature fluctuations in the cosmic microwave back- 
ground radiation arising from the density fluctuations re- 
sponsible for the structure formation in the universe. It 
was argued in ref. Q that it is possible to get linhm as 
large as few meters. 

There are many possible sources of density fluctuations 
in the early universe, inflation being one such possibility. 
Cosmic defects (strings, monopoles, textures) also lead 
to density fluctuations in the universe. There has been 
extensive study of density fluctuations generated by cos- 
mic strings from the point of view of structure forma- 
tion M . In this paper we study the influence of density 
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fluctuations produced by cosmic strings on the dynam- 
ics of the quark-hadron transition. Though recent mea- 
surements of temperature anisotropies in the microwave 
background by BOOMERANG, and MAXIMA experi- 
ments H at angular scales of i ~ 200 disfavor models of 
structure formation based exclusively on cosmic strings 
Still, due to many uncertainties in the scaling 
models of cosmic string network evolution one can not 
rule them out as candidates of sources of required den- 
sity fluctuations. Some feature of the evolving string net- 
work which is capable of introducing super horizon string 
correlations, may be able to make these models compati- 
ble with recent observations. Further, even with present 
models, it is not ruled out that cosmic strings may con- 
tribute to some part in the structure formation in the 
universe. Above all, cosmic strings generically arise in 
many Grand Unified Theory (GUT) models. If the GUT 
scale is somewhat lower than 10 16 GeV then the result- 
ing cosmic strings will not be relevant for structure for- 
mation, but they will still affect various stages of the 
evolution of the universe in important ways. (For a dis- 
cussion of these issues, see [0.) It is in this spirit that 
we undertake the study of the effects of cosmic strings on 
quark-hadron phase transition. 

We will show that very interesting effects can arise due 
to cosmic strings, such as generation of sheet like bary- 
onic inhomogeneities spreading across the horizon. Typ- 
ical separation between such sheets will be of the order 
of a km at quark-hadron transition stage. (For simplicity 
we do not consider here the effects of density fluctuations 
due to oscillating string loops. These will give rise to 
localized regions of baryon inhomogeneity. ) There is no 
other model known which can lead to such large values of 
separations between baryon inhomogeneity regions. Such 
inhomogeneities will have important effects on the nu- 
cleosynthesis. One could also reverse the argument and 
study the constraints on various cosmic string parame- 
ters which will arise from various elemental abundances 
from inhomogeneous nucleosynthesis in such a model. 

As we mentioned above, a re-consideration of the mod- 
els of IBBN is well motivated at the present stage due 
to recent measurements of the angular power spectrum 
of CMBR jij. These results imply a value of which 
is larger than that allowed by the conventional homoge- 
neous nucleosynthesis. In this context we mention that 
a recent analysis by Kurki-Suonio and Sihvola Q shows 
that the results of IBBN may be (marginally) compati- 
ble with this large value of Qj, as required by the recent 
CMBR measurements Q. However, as pointed out in 
rcf . Q , this agreement can be achieved only if regions of 
baryon inhomogeneity are separated by a distance scale 
of about 35-70 km (at T = 1 MeV), and if these re- 
gions have a more planar structure with high surface to 
volume ratio. It is intriguing that the baryon inhomo- 
geneities generated by cosmic strings in our model satisfy 
both these criterion. These are clearly planar structures. 
Further, distance between these sheets at quark-hadron 
transition being of order 1 km (or somewhat smaller) will 



translate to about 100 km distance scale at T = 1 MeV. 



II. COSMIC STRING WAKES IN A 
RELATIVISTIC FLUID 

We first briefly review the structure of density fluc- 
tuations produced by a cosmic string moving through a 
relativistic fluid. The space-time around a straight cos- 
mic string (along the z axis) is given by the following 
metric |L2|, 

ds 2 = dt 2 - dz 2 - dr 2 - (1 - AGp) 2 r 2 d<j) 2 , (1) 

where p, is the string tension. This metric describes a 
conical space, with a deficit angle of 8irGp. This met- 
ric can be put in the form of the Minkowski metric by 
defining angle <p' = (1 — AGp)cf>. However, now <f>' varies 
between and (1 — 4Gp)2ir, that is, a wedge of opening 
angle SirGp is removed from the Minkowski space, with 
the two boundaries of the wedge being identified. It is 
well known that in this space-time, two geodesies going 
along the opposite sides of the string, bend towards each 
other [ pH] , This results in binary images of distant ob- 
jects, and can lead to planar density fluctuations. These 
wakes arise as the string moves through the background 
medium, giving rise to velocity impulse for the particles 
in the direction of the surface swept by the moving string. 
For collisionless cold dark matter particles the resulting 
velocity impulse is Q, v mp i s — 4:irGpv st j s t (where v st 
is the transverse velocity of the string). This leads to a 
wedge like region of overdensity, with the wedge angle 
being of order of the deficit angle, i.e. 8irGp (~ 10~ D 
for GUT strings). The density fluctuation in the wake is 
of order one. Subsequent growth of this wake by gravi- 
tational instability has been analyzed in great detail in 
literature Q. 

The structure of this wake is easy to see for collisionless 
particles (whether non- relativistic, or relativistic). Each 
particle trajectory passing by the string bends by an an- 
gle of order AirGp towards the string. In the string rest 
frame, take the string to be at the origin, aligned along 
the z axis, such that the particles are moving along the 
—x axis. Then it is easy to see that particles coming 
from positive x axis in the upper/lower half plane will all 
be above/below the line making an angle ^fAnGp from 
the negative x axis. This implies that the particles will 
overlap in the wedge of angle 8irGp behind the string 
leading to a wake with density twice of the background 
density. One thus expects a wake with half angle W and 
an overdensity Sp/p where fl4| ], 

6 W ~ AirGp, ^ ~ 1. (2) 
P 

However, our interest is in the density fluctuations gen- 
erated by the strings at the quark-hadron transition in 
the early universe. At that stage, it is not proper to take 
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the matter as consisting of collisionless particles. A suit- 
able description of matter at that stage is in terms of a 
relativistic fluid which we will take to be an ideal fluid 
consisting of elementary particles in the QGP phase, and 
consisting of hadrons in the hadronic phase. Genera- 
tion of density fluctuations due to a cosmic string mov- 
ing through a relativistic fluid has been analyzed in the 
literature [p~5|— p~T|| . The study in ref . jL5| focussed on the 
properties of shock formed due to supersonic motion of 
the string through the fluid. In the weak shock approxi- 
mation, one finds a wake of overdensity behind the string. 
In this treatment one can not get very strong shocks with 
large overdensities. In refs. [ p3[|l^ ], a general relativistic 
treatment of the shock was given which is also applicable 
for ultra-relativistic string velocities. The treatment in 
ref. ]l7| is more complete in the sense that the equations 
of motion of a relativistic fluid are solved in the string 
space-time (Eq.(l)), and both subsonic and supersonic 
flows are analyzed. One finds that for supersonic flow, 
a shock develops behind the string, just as in the study 
of ref. (l^Jlfl ■ In the treatment of ref. jl^] one recovers 
the usual wake structure of overdensity (with the wake 
angle being of order Gp) as the string approaches ultra- 
relativistic velocities. Also the overdensity becomes of 
order one in this regime. For smaller string velocities, 
the results of M| and Jl(],[n| are roughly in agreement. 

In the following we will reproduce some of the relevant 
results from these papers. It is not expected that the 
string will move with ultra-relativistic velocities in the 
early universe. Various simulations have shown that rms 
velocity of string segments is about 0.6 jl8| for which 
the shock will be weak. In this situation, any of the ap- 
proaches of refs. p^-[l7|] can be followed. We will briefly 
discuss the properties of shock, along the lines of discus- 
sion in ref. |15|. We will also reproduce some relevant 
results from ref. [[t7| . These two sets of results will span 
the range of various parameters of interest for us. 

The basic idea in the study of ref. ]l5| is that the situ- 
ation of a string moving through the fluid is the same as 
the situation when the fluid flows past a wedge with the 
wedge angle being the same as the deficit angle associ- 
ated with the cosmic string. If the string velocity exceeds 
the sound speed in the fluid then one gets shock waves. 
The properties of the shock can be determined following 
the standard treatment of shocks |19|. 

Fig.l shows the two dimensional spatial slice, with the 
string being perpendicular to the plane and located at 
the point marked as O. Direction of fluid flow is shown 
by arrows. Dark shaded region denotes the region of 
the shock with overdensity. The fluid region in the front 
of the shock has background energy density. The wedge 
with white interior denotes the deficit space, with the up- 
per edge identified with the lower edge. Due to reflection 
symmetry of the problem, and due to the identification 
of upper and lower wedge boundaries, it is easy to see 
||15t that the direction of flow must be parallel to the 
upper and lower wedge boundaries respectively (forming 
the deficit region). The situation is thus indistinguishable 



from the situation of fluid flow past a wedge with wedge 
angle being the same as the deficit angle, i.e. 8wGp. We 
will denote the fluid variables in front of the shock and 
in the back of the shock by using subscripts 1 and 2 re- 
spectively. 




FIG. 1. 

Fluid flow past the cosmic string. Dark shaded re- 
gion denotes the region of the shock with overdensity. 
The fluid region in the front of the shock has background 
energy density. The wedge with white interior denotes 
the deficit space, with the upper edge identified with the 
lower edge. Arrows denote the direction of fluid flow. 

The energy momentum tensor of the fluid is taken to 
be of the perfect fluid form: 



= wu^u" + pr\ 



(3) 



Here w = p + p, p and p are energy density and pres- 
sure respectively, p, v = 0,1,2,3. Note that we use the 
flat metric rf v here as the problem of conical space-time 
is now mapped to that of fluid flow past a wedge in the 
Minkowski space-time. 

Certain continuity conditions must be satisfied at the 
surface of the shock. These arc [H , 



rpOn rpOn 

1 l 2 ' 



(4) 



Here superscript n denotes the direction normal to the 
shock surface and subscripts 1 and 2 denote pre-shock 
and post-shock regions. These conditions give (using 
Eq.(3)), 



7l„(Pl U ln +Pl) = llniP2vl n +P2)- 



(5) 
(6) 



Here 7i„ = 1/ yl — v\ n is the relativistic 7 factor cor- 
responding to velocity v\ n (similarly for 72^). In the rest 
frame of the shock surface, we have, 



Vln = Vfi 



AnGp), 



V2n = v 2 sin(6). (7) 



Vf is the velocity of the fluid in the cosmic string rest 
frame in region 1. Thus, Vf is the string velocity v string 
in the background plasma rest frame. Now, in addition 
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to the continuity conditions in Eqs.(5),(6), we also re- 
quire continuity of the component of velocity parallel to 
the shock. That is, 



V2cos(0) = vicos(0 + AirGp). 



(8) 



Plots of 8 (in degrees), and Sp/p, for a range of string 
velocities as obtained from Eq.(9) (solid curves) and 
Eqs.(lO) (dashed curves). Thick curves denote the val- 
ues corresponding to v s = l/y3, while thin curves cor- 
respond to v s = 0.1 



By writing the overdensity Ap = p2 ~ pi and over- 
pressure Ap = P2 — Pi, (and using Gp, << 1), one can 
determine various properties of the shock using above set 
of equations. We get, 



Sp 
Pi 



WnGpv'j 



sinO 



(9) 



For the first equation above we have used p\ = pi/3 
appropriate for relativistic ideal fluid. For the second 
equation we have used = where v s — is 
the isentropic sound speed. (To show that the shock is 
isentropic, one needs to consider additional continuity 
condition for particle number, see |jl5| , ^9| ). 

Eq.(9) contains all the relevant information for the 
overdense region behind the string. The regime of va- 
lidity of these equations is Ap/p, Ap/p << 1. For ultra- 
relativistic case, we use expressions from ref. Jl7[ . These 
are expressed in terms of fluid and sound four velocities, 
(in weak shock limit), 



Sp 16-!rGpuj(l 



Pi 



3u. t 



7 



sinO ~ — , 



(10) 



where Uf — vj/^Jl — vj and u s = v s / \/l — We see 

that in the non-relativistic regime Eq.(9) and Eq.(10) 
give roughly similar results. However, in the ultra- 
relativistic regime the shock properties obtained from 
these two sets of equations are very different (still within 
the regime of validity of the two treatments). The pa- 
rameters of relevance for us are 6 and Sp/p. In Fig. 2 we 
have given the plots of the values of these as obtained by 
the above two sets of equations. As we are interested in 
temperatures close to T c , in a close enough neighborhood 
of T c , the speed of sound becomes very small p0[ . Thus, 
in Fig. 2 we give plots for the value of speed of sound 
v s = l/\/3 (shown by thick curves), as well as v s = 0.1 
(thin curves) Q . 
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FIG. 2. 
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When string velocity is ultra-relativistic, then one can 
get strong overdensities (of order 1) and the angle of the 
wake approaches the deficit angle ~ 8nGp. For string ve- 
locities in the range of 0.6 - 0.9 the predictions of the two 
sets are roughly similar (within a factor of order one). We 
will use a sample value corresponding to string velocity 
of 0.9 for which we take, 



20°, Sp/p ~ 3 x 10" 



(11) 



These values correspond to those obtained from Eq.10 
for v s = 1/V3- As mentioned earlier, the overdensities 
expected of collisionless matter are of order 1, with the 
wake angle being of order of the deficit angle (as given in 
Eq.(2)). This is what one also gets from the treatment in 
]l7t for ultra-relativistic string velocities. Thus we will 
also consider this set of values. With this, we are now set 
to get back to our problem of understanding the effects 
of these density fluctuations on quark-hadron transition 
in the universe. 



III. EFFECT OF STRING WAKES ON 
QUARK-HADRON TRANSITION 

In the conventional picture of the quark-hadron tran- 
sition, the transition proceeds as follows |l| . As the uni- 
verse cools below the critical temperature T c of the tran- 
sition, hadronic bubbles of size larger than a certain crit- 
ical size can nucleate in the QGP background. These 
bubbles will then grow, coalesce, and eventually convert 
the QGP phase to the hadronic phase. The details of 
this transition lead to an interesting and complex pic- 
ture §]2C] HHH. Very close to T c the critical size of the 
bubbles is too large, and their nucleation rate too small, 
to be relevant for the transition. Universe must super- 
cool down to a temperature T sc when the nucleation rate 
becomes significant. The actual duration of supercooling 
depends on various parameters such as the values of sur- 
face tension a and the latent heat L. We take the values 
of these parameters as in ref. || (motivated by lattice 
simulations |J]), a = 0.015T C 3 and L ~ 3T C 4 . With these 
values, one can estimate the amount of supercooling to 
be ||| (we take T c = 150 MeV), 



AT SC = 1 - ^ 



10" 



(12) 



As the universe cools below T sc , bubbles keep getting 
nucleated and keep expanding. This nucleation process 
is very rapid and lasts only for a temperature range of 
AT n ~ 10~ 6 , for a time duration of order At n ~ 10 _5 £# 
(tn is the Hubble time) p3"P]. The latent heat released 
in the process of bubble expansion re-heats the universe. 
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Eventually, the universe is reheated enough so that no 
further nucleation can take place. Further conversion of 
the QGP phase to the hadronic phase happens only by 
the expansion of bubbles which have been already nu- 
cleated. Even this expansion is controlled by how the 
latent heat is dissipated away from the bubble walls. Es- 
sentially, the universe cools little bit, allowing bubbles to 
expand and release more latent heat. After the phase of 
rapid bubble nucleation, the universe enters into the slow 
combustion phase Jl| . As we will show below, the picture 
of this slow combustion phase is very different in our 
model. Here, in the wake region of the cosmic string, the 
conversion of the QGP phase to the hadronic phase pro- 
ceeds by moving macroscopic planar boundaries. First, 
let us recall the estimate of the duration of this slow 
combustion phase. 

Energy density, entropy density, and pressure 
(pq,s q ,p q ) in the QGP phase and (ph,Sh,Ph) in the 
hadronic phase are, 

Pq = g q aT 4 + B, s q = ^g q aT 3 , Vq = ^aT 4 - B, (13) 

p h = gnaT 4 , s h = ^g h aT 3 , p h = yaT 4 . (14) 

Here g q ~ 51 and gu — 17 are the degrees of freedom 
relevant for the two phases respectively (taking two mass- 
less quark flavors in the QGP phase, and counting other 
light particles) pl{ |. At the transition temperature, we 
must have p q = ph which relates T c and the bag constant 
B, 



B = ^a T c(9g ~ 9h)- 



(15) 



For the transition occurring at the temperature T ~ 
T c , the latent heat is, L = T c (s q — s/,) = AB. For the time 
evolution of the scale factor R of the universe, Einstein's 
equations give, 



R(t) 



8nGp(t) _ [J 
3 V B' 



(16) 



where A -1 gives the typical time scale at the quark- 
hadron transition. For T c = 150 MeV, A^ 1 ~ QApsec. 
Now, conservation of the energy-momentum tensor gives, 



R(t) 



3 Mt) 



dt 



dt 



{R(tf[p(t)+p(t)}}. 



(17) 



For the mixed phase, we write the average value of 
energy density as, p — fp q + (1 — ,f)ph, where / is the 
fraction of volume in the QGP phase. Now, during the 
transition, T and p are approximately constant. Thus 
Eq.(17) gives, 



R _ -[f(Sq - Sh) 

R 3[f(s q - s h ) + s h ] ' 



(18) 



From Eqs.(16),(18), we can write down the evolution 
equation for the fraction /, 



• /+ 0^T) [4/ + ^T ]1/2{3/(x " 1) + 3} = ' (19) 



where we have defined \ = ff" — 3. From this equation 
one can estimate the total duration of the quark-hadron 
transition Ai trnsn , that is the time in which the fraction 
/ decreases from initial value 1 (QGP phase) to final 
value (hadronic phase). We find that Attmsn — 14/x 
sec pj| . Compared to this the Hubble scale at the time 
when T = T c is tn — 50/i sec. 

Let us now see how density fluctuations in the shocks 
created by cosmic strings affect the dynamics of quark- 
hadron transition. We first take the parameters of the 
density fluctuations as given in Eq.(ll). Density fluctu- 
ation Sp/p ~ 3 x 10~ 5 translates to a temperature fluc- 
tuation of order AT wake = ST/T ~ 10~ 5 . We note that 
this temperature fluctuation is larger than AT n ~ 10~ 6 . 
This means that the nucleation of hadronic bubbles will 
get completed in the QGP region outside the wake of the 
string, while no nucleation can take place in the wake 
region during that stage. Since AT SC is of similar order 
as AT wa ke, one can conclude that the outside region will 
enter a slow combustion phase before any bubbles can nu- 
cleate in the region of overdensity in the wake. For this 
it is required that the overdensity in the wake should not 
decrease at the time scale of At n ~ 10 _5 tjy. The typical 
average thickness of the overdense region in the shock 
d s hk w iU be (for a wake extending across the horizon), 



* sh k 



sinO m — 5km. 



(20) 



A density fluctuation of this length scale (which is 
smaller than the horizon size) primarily propagates as 
a plane wave (though, since d s hk is not too small com- 
pared to r/i, the evolution of this overdense region will 
have two modes, one propagating like plane wave and 
the other growing as R 2 |2jJ). Typical time scale for the 
evolution of the overdensity in this region, t s hk , will be 
governed by the sound velocity v s which becomes very 
small close to the transition temperature dJ|^]. Taking 
v s ~ 0.1 we get, 



ishk 



^shk 



t H . 



(21) 



Thus the density (and hence temperature) evolution 
in the shock region happens in a time scale which is too 
large compared to At n . It is also larger (by a factor of 
about 4) than At tr nsn during which the quark phase is 
completely converted to the hadronic phase in the region 
outside the wake. We mention that in our picture, we 
consider the time when the universe has just started go- 
ing through the quark-hadron transition, and we focus 
on a region in which a wake of density fluctuation of size 
of order th has been created by the moving string. Es- 
sentially the region of study for us is the horizon volume 
from which the string is just exiting at the time when the 
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universe temperature T — T sc . The formation of most of 
the region of shock thus happens when the temperature 
is still large enough compared to T c so that the speed of 
sound is close to the value However, some portion 

of the wake will certainly form when the temperature 
is close enough to T c , that the relevant sound speed is 
small, say, v s ~ 0.1. At that stage the parameters of 
the shock will be determined by the plots shown by thin 
curves in Fig. 2. The extent of the overdense region will 
be governed by the time scale t s hk- Thus, if t s hk is much 
smaller than tn, then wake will not extend across the 
horizon. 

The precise time duration, Ati ag , by which the process 
of bubble nucleation in the shock region lags behind that 
in the region outside the wake is given by (|J2^] , 

At lag ~ ATw 3 a v f H c 10-H H , (22) 

with v s = 1/vo. Atiag will be larger if we take v s = 0.1. 
Ati a g is the extra time in which the temperature in the 
wake decreases to T sc compared to the time when the 
temperature drops to T — T sc in the region outside the 
wake. Since Ati ag is of same order as At n , we conclude 
that the region outside the wake enters the slow combus- 
tion phase before any bubbles can nucleate in the wake 
region. It is then reasonable to conclude that the la- 
tent heat released in the region outside the wake will 
suppress any bubble nucleation in the wake especially 
near the boundaries of the wake. If the heat propagates 
to the interior of the wake, then the bubble suppression 
may extend to the interior of the wake also, implying 
that there will simply be no bubbles in the entire wake 
region. In that case, hadronic bubbles which have been 
nucleated outside the wake will all coalesce and convert 
the entire outside region to the hadronic phase (with oc- 
casional QGP localized regions embedded in it). This 
hadronic region will be separated from the QGP region 
inside the wake by the interfaces at the boundaries of the 
wake, as shown in Fig. 3a. 

Further completion of the phase transition will happen 
when these interfaces move inward from the wake bound- 
aries. Clearly, these moving, macroscopic interfaces will 
trap most of the baryon number in the entire region of the 
wake (and some neighborhood) towards the inner region 
of the wake. Finally the interfaces will merge, complet- 
ing the phase transition, and leading to a sheet of very 
large baryon number density, extending across the hori- 
zon. Actual value of baryon density in these sheets will 
depend on what fraction of baryon number is trapped in 
the QGP phase by moving interfaces. It is possible that a 
large fraction of the total baryon number in the wake re- 
gion may get trapped in these sheets, resulting in sheets 
of extremely high baryon densities [gj . The baryon num- 
ber density on the sheet may vary by a factor of two, from 
one end of the sheet to the other end, since the shock re- 
gion has the form of a wedge (though, the portion of the 
shock region formed later may have different parameters 



as compared to the portion formed in the beginning due 
to different sound speeds). 

Second possibility is that the bubble nucleation is 
not entirely suppressed near the center of the shock re- 
gion. While the region outside the wake converts to the 
hadronic phase, same may happen at the center of the 
wake as well. In that case the hadronic phase will spread 
from inside the wake at the same time when the hadronic 
phase is moving in from outside the wake through the 
wake boundary. These two sets of interfaces will then 
lead to concentration of baryon number in two different 
sheet like regions, with the separation between the two 
sheets being of the order of a km or so. This type of 
evolution is shown in Fig. 3b. (Note, from Eq.(19) that 
the total time duration of the transition still remains the 
same, irrespective of the fact that now the transition is 
happening by moving planar boundaries, rather than by 
expansion of bubbles.) The final result is that, string 
induced shocks will lead to either a single sheet of high 
baryon concentration spreading across the horizon, or a 
pair of sheets again spreading across the horizon. 

Typical separation between such structures will be gov- 
erned by the number of long strings in a given horizon, 
which is expected to be about 15 (from numerical simu- 
lations (26|). Thus these sheets (or pairs of sheets) will 
be separated by a typical distance of order of a km which 
is far too large compared to any other relevant scale at 
that time which could dissipate such inhomogeneities || . 
These structures will then certainly survive until the time 
of nucleosynthesis and will affect the abundances of ele- 
ments K|. We mention here that we are not consider- 
ing the effect of density fluctuations produced by string 
loops. These will also lead to baryon number inhomo- 
geneities via the effects discussed here. However, these 
structures will be on a more localized scale. It is more 
complicated to calculate the effects of density fluctua- 
tions by oscillating loops (especially when time scales are 
of crucial importance). Still, a more complete investi- 
gation of the effects of cosmic strings on quark-hadron 
transition should include this contribution also. 




(a) (b) 



FIG. 3. 
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These figures show portions of the overdense shock re- 
gion. Q and H denote QGP and hadronic phases respec- 
tively and solid lines denote the interfaces separating the 
two phases, (a) Situation when no hadronic bubbles nu- 
cleate inside the overdense wake region, (b) If hadronic 
bubbles can nucleate in the inner region of the wake also 
then there will be two sets of interfaces trapping QGP 
regions in two separate sheets. Shaded regions, marked 
as B, denote sheet like baryon inhomogeneities expected 
in the two cases. 

We now discuss the situation if the shock behind the 
string is strong. In that situation one expects a wake of 
opening angle of about 10 -5 with overdensities of order 
1. These overdensities will lead to temperature fluctu- 
ations of order 1 behind the string. This temperature 
fluctuation is far too large compared to the temperature 
interval of bubble nucleation AT n , or even AT SC . In such 
a situation, even if the transition to the hadronic phase 
has been (just) completed, as string moves through the 
hot hadronic medium, the overdensity in the wake will 
force the temperature in the wake region to rise above 
T c , forcing that region to convert to the QGP phase. The 
evolution of this region will be complicated due to the 
tension of the interfaces at the boundaries of this region. 
(Otherwise an overdensity of wavelengths much smaller 
than the horizon simply propagates as plane waves). If 
we assume that the region primarily cools due to the uni- 
verse expansion, then a wake of QGP region extending 
across the horizon will form before the interfaces close 
in and convert this temporary QGP phase back to the 
hadronic phase. End result will again be a sheet of large 
baryon density spreading across the horizon. (In this 
case, the separation between the two sheets for the case 
of Fig. 3b will be of order few cm, as the average thickness 
of the wake is about 0.5 x 10 _5 x 15 km ~ 8 cm. Thus, 
by the time of nucleosynthesis, effectively only one sheet 
may survive. Also, due to smaller wake thickness, baryon 
density on these sheets in this case will be smaller.) If 
the region cools at a faster rate then the sheets will be 
of smaller length. 



to horizon size sheets, or pairs of sheets with large baryon 
number densities. These will have strong effects on the 
abundances of elements as calculated during the epoch 
of nucleosynthesis. This becomes important in view of 
the recent measurements of the angular power spectrum 
of CMBR |I| which imply a large value of fib which is 
not compatible with the calculations of conventional ho- 
mogeneous nucleosynthesis. Our model leads to a very 
specific structure of baryon inhomogeneities. These are 
planar structure (with large surface to volume ratio), and 
are separated by a distance scale of about 1 km at the 
quark-hadron transition (which translates to a distance 
scale of about 100 km when T = 1 MeV). Baryon inho- 
mogeneities of this nature seem to be what is required to 
make the results of IBBN (marginally) compatible with 
the large value of f2f, required by the recent CMBR mea- 
surements |Q. 

Observations of these abundances may also be used 
to constrain the cosmic string parameters as well as the 
models of cosmic string network evolution. Since the 
baryon number density in these sheets will be very large, 
it will be interesting to check whether such sheet like re- 
gions can survive up to the present stage. Finally we 
mention that the effects described here can lead to in- 
teresting consequences for electroweak (EW) phase tran- 
sition. For example, just after the EW transition, EW 
symmetry may get restored in the overdense wake region 
behind the cosmic string. For a strong shock, that re- 
gion may cool rapidly, breaking the EW symmetry spon- 
taneously and generating baryon number in the process 
(even for a 2nd order EW phase transition). We will 
discuss this in a future work. 
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IV. CONCLUSION 



We conclude by re-emphasizing that cosmic string de- 
fects are generic in GUT models, irrespective of their role 
in structure formation and it is important to study how 
they can affect the dynamics of the universe. Dynamics 
of phase transitions, especially for a first order transi- 
tion, can be extremely sensitive to temperature varia- 
tions. Density fluctuations, and hence temperature fluc- 
tuations, produced due to cosmic strings, even if they are 
of small magnitude, can affect the phase transition dy- 
namics in crucial ways. In this paper we have shown that 
the presence of cosmic strings can dramatically alter the 
dynamics of a first order quark-hadron transition leading 
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